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ABSTRACT 

Motivated both by considerations of the generation of large-scale astrophysical 
magnetic fields and by potential problems with mean magnetic field generation by 
turbulent convection, we investigate the mean electromotive force (emf ) resulting 
from the magnetic buoyancy instability of a rotating layer of stratified magnetic 
field, considering both unidirectional and sheared fields. We discuss why the 
traditional decomposition into a and j3 effects is inappropriate in this case, and 
that it is only consideration of the entire mean emf that is meaningful. By 
considering a weighted average of the unstable linear eigenmodes, and averaging 
over the horizontal plane, we obtain depth- dependent emfs. For the simplified 
case of isothermal, ideal MHD we are able to obtain an analytic expression for 
the emf; more generally the emf has to be determined numerically. We calculate 
how the emf depends on the various parameters of the problem, particularly the 
rotation rate and the latitude of the magnetic layer. 

Subject headings: Sun: interior — Sun: magnetic fields — instabilities — MHD 
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1. Introduction 

One of the most important theoretical problems in astrophysical fluid dynamics is to 
explain the generation of global scale magnetic fields, as observed in stars, accretion discs 
and galaxies. It is generally accepted that these are the result of some sort of hydromagnetic 
dynamo, in which the inductive motions of a highly electrically conducting plasma maintain 
magnetic fields against their tendency otherwise to decay. However, a convincing theoretical 
explanation of how such large-scale magnetic fields can be generated — i.e. fields with a 
significant component on scales much larger than that of the plasma motions responsible 
for their generation — remains elusive. 



The traditional theoretical approach to explaining the generation of 



fields is via mean field electrodynamics (see, for example, 



Moffatt 



1978 



arge-scale magneti c 



Krause fc Radler 



1980|) , an elegant theory of magnetohydrodynamic (MHD) turbulence in which the evolution 
of the mean (large-scale) field is governed by a mean induction equation of the form 

—1 = V x (U x B ) + V x £ + V V 2 B , (1) 

where B Q represents the mean magnetic field, Uq the mean velocity, £ the mean 
electromotive force (emf) and rj the magnetic diffusivity. The term describing the mean 
emf, which is the distinctive feature of equation (CTJ in comparison with the unaveraged 
induction equation, is defined by 

£ = (u x b), (2) 

where u and b represent the (small-scale) fluctuating velocity and magnetic fields, and 
angle brackets denote a spatial average over intermediate scales. The closure of equation (CTJ 
is usually brought about by postulating an expansion of £ in terms of Bq and its spatial 
derivatives, 

Si = aijB j + fiijk + • • • , (3) 
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where a. and j3 are pseudo-tensors; iHughes fc Proctor! (120101 ) discuss a more general 



expansion procedure involving also temporal derivatives of the mean field. In the kinematic 
regime, in which the field is assumed to exert no back-reaction on the flow, the components 
aij and depend solely on the properties of the velocity field and on the magnetic 
diffusivity. The symmetric part of the a tensor (the so-called 'a-effect') leads to field 
amplification, and can be non-zero only in flows that lack reflectional symmetry, such as 
helical flows; in its simplest isotropic form, in which = /Se^fe, the scalar /3 can be 
identified as a turb ulent diffusivity, althou gh in general it has a much more complicated 



interpretation (see 



Krause fc Radler 



19801 ) . Often, astrophysical magnetic fields are 
modelled by assuming plausible (although somewhat arbitrary) functional forms for a, 
(3 and Uq, including the parameterisation of nonlinear effects. Via this procedure, it is 
possible to reproduce many of the features of observed magnetic fields, such as the solar 
cycle. 

However, notwithstanding the ability of mean field models to be able to mimic cosmic 
magnetic fields, there remain some potentially serious difficulties regarding the application 
of the standard formulation of mean field electrodynamics to turbulent flows at high 
magnetic Reynolds numbers, the case of astrophysical relevance. Of particular significance 
for the type of model we shall consider in this paper are those that may be encountered 
in the kinematic regime. When the magnetic Reynolds number Rm is small, everything 
is fine; magnetic fluctuations arise only from the interaction between the mean field and 
the flow, b is then linearly related to the mean field Bq, the expansion ([3]) is well-founded, 
and equation (TjQ) provides the correct description of large-scale magnetic fields. However, 
under the astrophysical conditions of high Rm turbulence, there is the likelihood that 
small-scale magnetic fields can be generated ind ependently of any large-sc ale mean field 



small-scale dynamo action (see, for example, 



Childress &: Gilbert 



19951 ). In this case 



there will be a competition between two distinct dynamo modes. The physical principles 
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underlying large- and small-scale dynamo action are very different: the former requires 
a lack of reflectional symmetry in the velocity (as manifested, for example, by the flow 
helicity) and has a dynamo growth rate that scales with the (small) wave number of the 
large-scale field; the latter has a growth rate that can be related to the chaotic properties 



of th e flow, characterised by Lyapunov exponents and cancellation exponents (]Du Ott 



19931 ) . The nature of t he fie ld th at would be observed in thi s case remains a controversial 



issue. 



Boldyrev et al. 



fl2005f ) and 



Cattaneo &: Hughes! (120091 ) have argued that at high 



Rm, when the turbulence is rotationally influenced but not rotationally dominant, and 
in the absence of a large-scale velocity shear, the small-scale mode will dominate; the 
observed fields will then bear no relation to those predicted from any calculation of cx 



and (3 and subsequent solution of equation jTJ. Numerical evidence 



has 



) een provided by simulation s of rotating convective turbul ence ( 



2006 



Hughes fc Cattaneo 



or this point of view 



Cattaneo fc Hughes 



20081 ) and forced helical turbulence ( ILivermore et al. 



20071 ). 



Conversely, if convection is strongly constrai ned by rotation then l arge-s cale field generation 
may be possible; this was first examined by 



numer i cal investigati o ns by , for example 



(120041 ) 



Kapyla et al. 



by 


Childress & Soward 


1972 


), with subsequent 


Jones & Roberts 


( 


2000) 


. Stellmach & Hansen 



( 120091 ). Whether the mean field idea holds at high Rm in the presence 



of a large-scale velocity shear remains a fascinating, but unresolved issue. Numerical 



simulations at moderate value s of Rm (see, for example, 



Yousef et al. 



2008 



2008 



Hughes &: Proctor 



Kapyla et al. 



20091 ) have shown the beneficial role of a large-scale shear flow for 



large-scale fi eld generation, but were n ot definitive in pinning down the underlying physical 



mechanism. iProctor fc Hughed (120111 ) have recently addressed this problem systematically 
by considering the nature of the magnetic fields generated by spatially-filtered flows driven 
by rotating convection. Intriguingly, the model offers support both to the standard mean 
field au picture as well as to the idea that large-scale field generation may come from 
large-scale flow components (a 'small scale' dynamo on the large scales). Further work is 
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taking place to understand this extremely important issue. 

Thus, given possible difficulties with the traditional kinematic a-effect picture, together 
with the intrinsic interest in alternative dynamo mechanisms, we explore in this paper the 
conceptually very different idea of an emf resulting from an instability of the magnetic field 
itself. Unlike the standard picture of a weak field being amplified by a velocity field that 
exists in the absence of a magnetic field — e.g. turbulence driven by thermal convection 

- here a pre-existing magnetic field is absolutely pivotal to the generation of an emf that 
can subsequently lead to field amplification. Such a means of generating an emf has been 



well-stud ied in p 



example, 



Capello 



asma confinement devices, particularly reversed field pinches (see, for 



2004]). It should though be stressed, since we are envisaging a dynamo 
process, that the idea of a magnetic field-driven emf does not of course provide 'something 
for nothing'. In reversed field pinches a strong toroidal field results from external coils. In 
any proposed astrophysical dynamo, it is envisaged that toroidal field is wound up from 
poloidal field by the differential rotation of the body, and that an instability of this field 
can then give rise to an emf that will close the dynamo loop by regenerating the poloidal 
field; the underlying energy source is thus the differential rotation, which itself may be a 
consequence of rotating thermal convection. 

Our work is particularly motivated by considerations of a possible dynamo located in 
the solar tachocline, with the idea that the strong velocity shear could generate toroidal field 
from a weak poloidal component, and in which the poloidal field could be maintained via 
an emf driven by magnetic buoyancy instability, the mechanism believed to be responsible 
for the escape of magnetic flux from the deep interior. The basic idea of magnetic buoyancy 
instability is that a horizontal magnetic field that decreases with height will 'puff up' the 
atmosphere, and that if the field gradient is sufficiently strong then instability can ensue. 
In its details though it can be quite a subtle instability, with 3D modes preferred, despite 
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their having to overcome the stabilising effects of magnetic tension, and with further 
complications arising from the roles of rotation and diffusion. A recent review of magnetic 
buoyancy instabilities can be found in iHughesi (120071 ) . 



At the onset of instability, the fastest growing linear eigenfunction will dominate, and 
hence, at least in the early stages of the instability, the ensuing emf will be driven by this 
linear mode. In this paper we therefore investigate the nature of the emf that results from 
the rotationally influenced magnetic buoyancy instability. We consider the instability in 
fairly general terms, retaining all possible modes of instability, and considering two types of 
horizontal equilibrium magnetic fields: unidirectional fields of the form B = (B x (z), 0, 0), 
and what in plasma physics are often referred to as 'sheare d' fields, takin g the form 



3 = ( B T (z) , B y (z), 0). We therefore build upon the work of 



(1984. 



Moffatt 



fll978h and 



Schmitt 



20031 ) who considered a similar 



instabilities of a unidirectional field. 



problem, but restricted attention to magnetostrophic 



Thelenl (120001 ) also calculated the emf resulting 



from a magnetic buoyancy instability, but considered the somewhat different case of 
Rayleigh- Taylor type instabilities driven by the density discontinuity at the upper interface 
of a layer of field embedded in a non-magnetic atmosphere. An ostensibly similar problem, 
though actually o ne wi th rather different physical foundations, has been considered by 



Ferriz-Mas et al. 



(120041 ) . who calculated the emf arising from the instability of isolated 



magnetic flux tubes; the differences between the instabilities of a continuo usly stra t ified 



field and those of isolated flux tubes, which are significant, a re dis cussed in 



Cline et al. 



Hughesl (120071 ) 



(120031 ) have shown, via fully 



In an interesting model with a different set-up, 
nonlinear numerical simulations, how dynamo action may be achieved by the combined 
action of a velocity shear and magnetic buoyancy; for their model, in contrast to that 
discussed in this paper, the flow and the instability are on the same scale, so a mean field 
picture is not appropriate. 
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The notion of an instability-driven dynamo has also been discussed with regard to 
other instability mechanisms, although not in terms of generating a mean electromotive 
force that can close the dynamo loop. The primary candidate for causing fluid turbulence 
in accretion disks, necessary for angular momentum transfer, is the magnetorotational 
instability (MRI), which, via only a weak magnetic field, c an lead to the instability of a 



differentially rotating disk that is hydrodynamically stable ( jBalbus fc Hawle 



199 lh . By a 



'bootstrapping' mechanism, the resulting three-dimensional velocity field can then act to 
amplify the magnetic fie ld that, via its instability, had created the velocity field originally 



( IBalbus & Hawley 



19981 ) . In terms of magnetic field generation, the role of the MRI is thus 
to transform the velocity field from a state of purely differential rotation, which cannot 
act as a dynamo, to a three-dimensional turbulent state, which can. However, exactly how 
this mechanism may work remains poorly understood. Many early numerical simulations 
produced healthy-looking dynamos, but often these employ ed only numerical dissipation. 



Fromang &: Papaloizou 



The dangers of such an approach have been highlighted by 
and it has become clear that if an understanding is to be obtained computationally 
high-resolution simulations with explicit dissipation will be needed. 



(I2007f ). 



Spruit 



(120021 ) has discussed, in general terms, how a dynamo might operate in 
convectively stable regions of stellar interiors, driven by differential rotation and a Tayler 
instability of the magnetic field (rather than magnetic buoyancy, as discussed here). Order 
of magnitude estimates are given of the strength of the magneti c field that may resu lt from 
such a mechanism, should it work. However, as pointed out by 



Mathis et al. 



( 120081 ). the 



closing of the dynamo loop in Spruit's model comes directly from the small-scale poloidal 
field, and not, as it should, from the large-scale emf resulting from the instability. This is an 
important point, since an instability per se is not enough, as can be seen by consideration of 
tw o- dimensional insta bilities; it is the mean emf that is crucial. The numerical simulations 
of iMathis et al.l (120081 ) reproduce the Tayler instability, but, in the parameter regimes they 
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explored, they do not find dynamo action. 



The underlying principles of instability-driven dynamos are very different from those 
of 'traditional' dynamo mechanisms. The latter are characterised by the exponential 
amplification of a weak seed magnetic field, which grows until the nonlinear influence 
of the Lorentz force back on the flow allows equilibration to a statistically steady state. 
The evolution of the field can be split into distinct kinematic and dynamic regimes; for 
the former, it makes sense to discuss properties of the velocity field itself that may be 
beneficial for magnetic field generation. Furthermore, for kinematic mean field dynamos 
one can decompose the mean electromotive force into distinct components proportional to 
the large-scale field and to its gradient, as suggested by equation d3J), although the caveats 
discussed above regarding the difficulties of this approach at high values of Rm need to be 
borne in mind. For instability-driven dynamos, the whole picture is very different since 
the perturbations to the flow and magnetic field needed to drive the emf come about only 
through an initial instability of the field itself. There is therefore no concept of a strictly 
kinematic regime in which the flow can be described independently of the magnetic field. 
Furthermore, although, as we shall show below, it is possible to describe such dynamos 
within the mean field framework, it is in general not meaningful to speak of a and /3-effects. 
The instability, and hence the resulting emf, will typically have a complicated nonlinear 
dependence on the strength and gradients of the magnetic field, rendering meaningless a 
decomposition such as ([3]). An exception to thi s might be the maKn etoro tational i nstab ility, 



which is a weak field instability. In this case, 



Brandenburg! (120051 ) and 



Gressej (120 loh 



have argued that a mean fiel d expans i on of the form (jBJ) is possible, although it should be 



stressed, as acknowledged by 



L0) 



Gressell (120101 ). that a more symmetric treatment that treats 



the momentum a n d indu ction equations on an equal footing, an approach pioneered by 



Courvoisier et al. 



(l2010al lbl). is required. 
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Obtaining a mathematical description of an instability-driven dynamo is less 
straightforward than for a standard, flow-driven dynamo, where, at least in the kinematic 
regime, the problem is tackled either via the induction equation or the mean induction 
equation. For an instability-driven dynamo, the starting point is an analysis of the pertinent 
instability. The difficulty then arises in describing how this instability feeds back into the 
dynamo process. The initial phases of the instability will be characterised by exponential 
growth, and it therefore makes sense as a first step to consider the nature of the emf (which 
will also be growing exponentially) resulting from the most unstable mode of instability and 
its dependence on the various parameters of the problem. The ultimate aim must be to try 
to incorporate, self-consistently, the emf arising from the instability, in both its linear and 
nonlinear phases, with the generation of toroidal field from differential rotation. 

Our aim in this paper is to evaluate the nature of the emf arising from magnetic 
buoyancy instability of a rotating, stratified layer of gas containing a horizontal, depth- 
dependent, but not necessarily unidirectional, magnetic field. In §[2] we describe the 
formulation of the instability problem and explain the prescription we use for calculating 
the resulting emf. In §[3] we consider the instability of a unidirectional field for the simplified 
case of isothermal, ideal MHD; this allows an analytical solution for the emf, and thus 
highlights its complex dependence on the magnetic field and its gradient. This also allows 
us to discuss the nature of the emf in the case of perfect electrical conductivity, an issue that 



19781 ). SectionS 



has given rise to problems in the standard mean field approach (iMoffattl 
considers the more general instability analysis, with all the diffusivities restored, which 
has to be performed numerically. In Section we extend our analysis to consider the emf 
resulting from the instability of a sheared magnetic field. Our conclusions and a discussion 
of how the work may be extended are contained in §[6j 
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2. Formulation of the Instability Problem 

We consider a plane layer of fluid of depth d located at colatitude 9, rotating with 
angular speed Q about the 6 = axis. We employ Cartesian coordinates, with x, y and z 
facing east, south and inwards respectively. The top of the layer is assumed to be the z = 
plane. The angular velocity is given by Q = Q(0, — sin#, — cos 8) = Qk. For simplicity, the 
fluid is assumed to be isothermal, with infinite thermal diffusivity. 

In standard notation, the governing equations of compressible MHD in a rotating frame 
may be expressed as 

+ 2p(n x u) = -Vp + —(V x B) x B + pge z + y. \V 2 u + ^V(V ■ u) j , (4) 
f) R 

V x (u x B) + 7]V 2 B, (5) 



Ot 



^ + P V ■ u = 0, (6) 
V • B = 0, (7) 

V = Pc\ (8) 

where c = (IZT) 1 ^ 2 is the isothermal sound speed, r\ is the magnetic diffusivity (assumed 
constant), \i is the viscosity (also assumed constant) and po is the magnetic permeability. 

We consider the instability of stationary equilibrium states, with uniform horizontal 
magnetic fields, linearly dependent on depth. In the stability analyses of §§EJH] the 
equilibrium field is assumed to be of the form B = B Q (1 + (z/d)e x , which we regard 
as the toroidal field. However, although it is thought that the toroidal field in the solar 
interior is dominant, it is of interest to consider the role of an additional weak poloidal 
ingredient; in §|5]we therefore extend our analysis to include equilibrium fields of the form 
B = B Q (1 + Cz/d)e x + B X {1 + &/d)e y . 
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Using equations (J4j) — (JH]) , and eliminating the pressure, linear perturbations to the 
basic state are governed by the following equations, here expressed in dimensionless form: 

p— + 2Qp (k x u) = - Vp + A 2 ((b • V)B + (B ■ V)6 - V(B ■ b)) + x pe z 
at 

+ z>(v 2 u+iv(V-u)), (9) 

- = Vx(tixB)+j)V 2 6, (10) 

^ + V-(pn) = 0, (11) 

V-b = 0, (12) 

where barred variables represent equilibrium values and unbarred variables represent 
perturbations. The equations are scaled using velocity c (the isothermal sound 
speed), length d (the layer depth), density po = p(z = 0) and magnetic field strength 
B xQ — B x (z — 0). The dimensionless parameters in equations - (fl2j) are defined by 

~ Qd 2 B 2 /p oPo gd fj, rj 

fi = — , A = 2 ' X=-2> v = j j r l = —>- ( 13 ) 

c c z p cd cd 

From now on we shall work only with dimensionless parameters, so, for simplicity, we shall 
drop the tildes. 

The equilibrium density distribution is calculated from the z-component of equation (J4]), 
after making use of (jSJ); in dimensionless form this may be written as 

Tz ~ x " = -~Tz B ■ (14) 

With the dimensionless equilibrium magnetic field written as B = (1 + (z)e x + A(l +£z)e y , 
the solution of (CHI) is given by 

+-(c 2 +a^ as) 
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where A = Bi/B (the ratio of the strengths of the fields in the y- and x-directions at 
z = 0). 

Since the equilibrium state is dependent only on z, the perturbations take the form 
^(a(z)(e^ kx+l ^ +st )), where a(z) is an amplitude function (which may be complex). The 
perturbation equations then become 

(4 \ v / klv\ A 2 A 2 

-k 2 + I 2 ) —u — ( 2f2 cos 9 H I v + 2VL sin 9 w + il — B v b x — ik — B v b v 
3 J i> V : '> l>) P P 

A 2 — i ik ikv . v .. 
+ —B x b z --p + —-w' + -u", 16 
P P 3 p p 

(fcl u \ / A \ u A 2 A 2 
2Q cos 9 — — — ) u — ( k 2 + -Z 2 ] — f — il—B x b x + ik—B x b v 
3 pj \ 3 J p p p 

A 2 -^'! H ilv . v .. 
+ —B y b z --p+--w'+-v", 17 
P p 3p p 



v A 2 — i A 2 —, A 2 — y 
sm) = -2Q sin 6> it - (A; 2 + / 2 ) -w - —B r b x - —Bh v + — (ikB x + UB V ) b z + -p 

p p p P P 

iku, ilv ft A 2 — , A 2 - 4i/_„ 



3p 3p p p p y 3 p 

sb x = ilByU - UB x v - B' x w - (A; 2 + l 2 )r]b x - B x w' + 776^, (19) 

sb y = -ikB y u + ikB x v - B~' y w - {k 2 + l 2 )i]b y - B y w' + 77&J, (20) 

s6 z = ^fcSj; + ZS 2/ )w - (A; 2 + l 2 )r]b z + r)b", (21) 

sp = —ikpu — ilpv — p'w — pw', (22) 

where B = (B x , B y , 0), and where a prime denotes differentiation with respect to z. 

An eigenvalue problem for the growth rate s, which may be complex, is formulated by 
dividing the layer between z = and z = 1 into 7N intervals (typically we take iV = 600), 
and using fourth order finite differences to approximate derivatives with respect to z. 
This allows the problem to be written in the form Ax. = sx, where A is a banded matrix 
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and where the eigenvector x takes the form (u , b , p , u\, &i, p%, . . . , u^, b^, p^). The 
eigenvalue problem is solved using the method of inverse iteration, which converges rapidly 
to the eigenvalue closest to an initial guess. 

The electromotive force for a particular mode of instability is determined by an average 
over the horizontal plane and, if the mode has a non-zero frequency, also over the period of 
oscillation. This lea d s to a ^-dependent mean emf, analogous in form to that considered by 



Childress fc Sowardl (119721 ) in their study of mean field generation in a rotating convective 
layer. The components of £ = {u x b) are given by 

£ x = l/2(v r b zr + Vib zi - w r b yr - Wib yi ), (23) 
E y = l/2(w r b xr + Wib xi - u r b zr - Uib zi ), (24) 
£ z = l/2(u r b yr + Uiy Z i - v r b xr - Vib xi ), (25) 

where u(z) = u r + iui, etc. As in any theory in which nonlinear terms are constructed 
from products of terms determined from a linear stability analysis, one has to deal with the 
issue of the amplitude of the terms, which, of course, is in some sense arbitrary. Here we 
choose to normalise the eigenfunctions with the total energy (kinetic + magnetic) of the 
perturbation. Since it is the curl of S that is of dynamical significance, rather than £ itself, 
then only £ x and £ y are of importance. 

A crucial point to recognise is that if there is more than one mode of in stability then 



1984 



20031 ). If there is a 



these will have different, and possibly competing, emfs ( ISchmittl 
unique mode of maximum growth rate then clearly the emf from this mode will dominate. 
If, however, there are two (or more) distinct modes of (equal) maximum growth rate 
then the emfs from each mode need to be considered. For example, in the non-rotating, 
unidirectional field case (i.e. when B y = and = 0), two modes with wave vectors (k,l) 
and (—A:, I) (where k, I > 0) have the same growth rate, but produce different emfs, with 
equal £ y but equal and opposite £ x . Thus, overall, there is no contribution to £ x , as might 
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be expected from standard mean field considerations for flows lacking reflectional symmetry 
(though it is important to remember that here we are dealing with the entire emf, not just 
the a-effect). Similarly, symmetry considerations suggest that, in the presence of rotation, 
£ x should vanish at the equator. 

Thus, in order to ensure that the emf obeys the symmetries of the problem, and also 
that it varies smoothly with parameters (for example, that 8 X = when Q = or at the 
equator), it is necessary to apply a weighting procedure between competing modes. Our 
prescription for the unidirectional field case is as follows. We determine the two modes of 
maximum growth rate with / > 0, one with k > and one with k < 0. If the two modes 
have growth rates Si and s 2 , with si > s 2 , and emfs £\ and £%, respectively, then we define 
the resulting emf by 

£ = £ 1 + e S2/si £ 2 . (26) 

Such a prescription allows a smooth transition as 9 passes through tt/2, for example. At 
the equator, modes with equal k, but oppositely-signed I have equal growth rates (si = s 2 ) 
and £ x and £ z of equal magnitude, but opposite sign. Away from the equator there is a 
definitely preferred mode. Expression ( 126]) incorporates these features. We employ a similar 
procedure for the instability of a sheared magnetic field, discussed in more detail in §[5j 

Finally in this section, we comment on the general form of the e mf that may be 



expec ted simply from considerations of the geometry of the problem. 



Krause fc Radler 



( 119801 ) derive an expression for £ in the case where there are three preferred directions, 
intended to represent the angular velocity f2, the mean velocity field U (which is zero in the 
present formulation) and the direction of gravitational attraction g, with the only unknowns 
being seve ral scalar quantities. As suming a linear dependence on each of the listed vector 



quantities, 



Krause fc Radlerl (jl980l . equation (5.45)), derive the following expression for the 
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mean emf: 

£ = -/3V xB-^gxB- i'U x B - f3 3 (O • V) B - /3 2 V (ft ■ U) - a[ (g ■ fl)B 

-a' 2 (g-B)fl-a' 3 (n-B)g. (27) 

Applying this to the formulation described above gives 
£ -- 



B (^f3 3 Q cos fl£ + a[Qg cos M + ^ j J + B x + jg (l + ^ j J ) e, 
+ cos fli + o/^ cos 0^1 + ^ j J - £? ^ + (l + ^ J J J e, (28) 



a 3 £l sin ^i?! ( 1 + -7 ] (?e 2 



d 

where a[, a 3 , (3, (3 3 and 7 are undetermined scalars. 

Owing to the influence of boundary conditions, we do not expect the emfs to have the 
exact form predicted by the above equations; we do however expect that the variation of 
the emf with the defining parameters of the problem will be similar. This can be measured 
by evaluating, for example, a root mean square emf defined by £ = (J £ 2 dz) . The above 
equations can be used to predict the dependence of £i on a selection of parameters, and 
these results compared with those obtained by solving the full eigenvalue problem. 



3. Ideal MHD Analysis of the Electromotive Force 



3.1. Calculation of the emf 



The first calculation o 
instabili t y was by 



Moffatt 



: the electromotive force resulting from magnetic buoyancy 



(119781 ). whose calculation was based on the stability analysis of 
Gilmanl (jl970[ ). with a small number of further approximations. In the analysis below, a 
si milar me t hod is followed, although no additional approximations are made beyond those 
of I Gilmanl ( I1970l ). We consider a rotating, isothermal layer of fluid, with sound speed 
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c, equilibrium magnetic fielc 
effects neglected. Following 



Be x and grav ity g = ge z , and with viscous and resistive 



Gilmanl (119701 ) . we shall consider only the component of 
rotation perpendicular to gravity — i.e. all terms involving cos 9 are neglected. As in §|2} 
perturbations are taken to be proportional to exp(i(kx + ly) + st)). The (dimensional) 
equations governing the perturbation amplitudes as functions of z (cf. equations ( JIBjl - 
( 1221) ) can then be written as 

(29) 
(30) 

' x + ikBb z + pg, (31) 

(32) 
(33) 
(34) 
(35) 
(36) 



spu = 


2VL sin 9pw — ikp + b z B , 


spv = 


—il(p + b x B) + ikbyB, 


spw = 


— 2f2 sin 9pu — p' — B b x — 


sb x = 


—ilBv — B w — Bw', 


sb y = 


ikBv, 


sb z = 


ikBw, 


sp = 


—ikpu — ilpv — pw — pw' 


V = 


pc 2 , 



The emf in the x-direction is given by £ x = l/23ft(t>6* — wb*). Using the above equations 
to evaluate this gives 



Sr~— I \U)\ $t(Q) + ' 



I Is 2 



2 dz 



where 



Q 



s 2 (c 2 (\n p)' + A 2 (ln B)') + i2kc 2 s£l sin 9 



(37) 



(3J 



s 2 (A 2 + c 2 ) + k 2 A 2 c 2 

s r = 3?(s) is the growth rate, and whe re we have used the fact that the most unstable mode 
has small k and large I ( )Gilmanlll970l ). The eigenvalue s depends on each of the variables 
in the basic state formulation, including the large-scale magnetic field B. Thus, even in 
this extremely simplified case, the emf cannot be meaningfully expressed in the form of 



-18 - 



expression ([3]). It is similarly unreasonable to expect that the emfs calculated for the more 
complex problems considered in §§IUE]can be decomposed in this way. 



A similar expression can be derived for S y = l/29?(u>&* — ub*). As shown by 



Gilman 



( 119701 ) . v and b y scale as 1 /I, leading to S x and £ z , but not £ y , scaling as k/l. Thus S y 
dominates in magnitude in this example, a feature that we may expect also to hold for the 
full problem incorporating diffusion. 



3.2. A well-defined emf in the absence of magnetic diffusion 

It is of interest, at this stage, to point out an important mathematical difference, for 
the case of a perfectly conducting fluid (rj = 0), between the standard kinematic formulation 
of the mean emf and that resulting from an MHD instability. In the traditional approach, 
a formal solution of the ideal (diffusionless) induction equation in terms of the Cauchy 
solution leads to expressions for and i n terms of int egrals along fluid particle 



pa ths of ay e raged quadratic Lagrangian quantities (jMoffatt 



by 



Moffatt 



3). 



197J. However, as pointed out 



(119741 ). there is no guarantee that these integrals will converge as t — » oo, and 
thus the status of the results is unclear. By contrast, for the instability- driven emf described 
above — and, we believe, for instability-driven emfs in general — there is nothing singular 
about the behaviour when 7] = 0. Expression (1371) . calculated with rj = 0, is well-defined 
and, were we to incorporate weak magnetic diffusion, the resulting emf would simply be a 
regular perturbation to the result (!3~Tj) . 

The key difference between the two systems lies in the spectra of the mathematical 
operators describing the two linear problems. For any non-zero value of rj, however small, 
the induction operator possesses a discrete spectrum, with well-defined eigenf unctions; 
however, for rj = the discrete spectrum disappears, together with the eigenf unctions. 
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Indeed, this singul arity in the induction ope rator lies at the heart of fast dynamo theory 



(see, for example, 



Childress fe Gilbert 



1995). By contrast, the linear operator for the 
stability problem has a discrete spectrum of unstable modes and associated eigenfunctions, 
leading to an unambiguous determination of the electromotive force. 



4. Unidirectional Magnetic Field 

In this section we consider the emf resulting from the instability of a unidirectional 
mean magnetic field of the form B = B Q (1 + (z/d)e x , with non-zero values of the viscosity 
and magnetic diffusivity. The instability evolves according to equations (jT6"l) - (l22l (with 
By = 0), and the mean emf is calculated via horizontal averaging, given by ( I23l) - (l25l) . and 
the weighting procedure fl26|) . 

Figure [T] shows the three components of £ plotted as functions of z, for a test case at 
colatitude 9 = 75°, representing the region in which most magnetic activity is observed on 
the solar surface. For this parameter set, the modes with largest growth rates have wave 
numbers (k,l) = (0.70367, 14.812) (with complex growth rate s = 0.082626 - 0.067823i) 
for k > 0, and (k,l) = (-0.70334, 14.812) (with s = 0.082602 + 0.067830i) for k < 0. The 
emfs in the figure are calculated from the mode with largest growth rate; when there is a 
unique mode of maximum growth rate, as here, the differences between the weighted and 
unweighted emfs are small. 

It is of interest to note that the ^-component of the emf (i.e. perpendicular to the 
imposed magnetic field) has the largest magnitude. The calculations in §[3] indicate that 
in the absence of viscosity and magnetic diffusivity, the ratio S x /S y is 0(k/l), which is 
vanishingly small for the preferred mode of I — > oo. Although now all diffusivities are 
restored, the physical mechanism underlying the instability is unchanged, so we may still 
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Fig. 1. — £ x , £ y and £ z versus depth, with 6 = 75°, Q = 0.1, A 2 = 0.1, ( = 1, v = f] = 
5 x 10" 5 , x = 2. 

expect the ratio of the emfs to be 0(k/l), and hence small. It can also be seen that the emf 
is localis ed towards the top of the layer. Consideration of the simplest instability criterion 



given by iGilmanl (119701 ). showing that a layer is unstable if the magnetic field strength 
decreases with height, yields for the current problem the condition £/(l + (z) > 0. The 
left-hand side of this inequality is clearly largest at the top of the layer, leading to the 
eigenf unctions (and hence the emfs) being peaked in this region. Figure [2] shows the x- and 
y-components of V x £; these are the dynamically significant elements in the evolution of 
the mean field, as explained in §[2j 




(Vx(?) s o 




Fig. 2. — The x- and y-components of V x £ for the same parameter values as in Fig. [TJ 



It is also of interest to examine how the emf varies with the parameters of the problem. 
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In order to do this and to be able to display the results in a meaningful way, we must 
first introduce an overall (^-independent) measure of the emf; with this in mind, we define 
£ where £{ — ij £\dz\ . Figures [3] - [6] show the variation of £ with fi, 9, ( and B 
respectively. Throughout all simulations, the dimensionless magnetic diffusivity, viscosity 
and gravitational field strength were kept fixed. The first two plots in each set of graphs 
show £ x and £ y as functions of the parameter under investigation. The third shows the 
growth rates of the fastest growing modes, optimised over k and I, in the k, I > (solid 
line) and k < 0, I > (dashed line) directions; as can be seen, these curves are essentially 
identical. Each graph is constructed using the weighting procedure defined by equation (T2T)]) . 
As discussed in §[2j this provides a sensible description of those cases for which two waves 
travelling in opposite directions produce cancelling emfs (such as at 9 = n/2). It can be 
seen that this feature is reflected in our results. 

Figure [3] shows the variation of £ x , £ y and growth rate s r with angular velocity Q. The 
symmetry considerations leading to equation f )28|) suggest that £ x oc Q and £ y = constant. 
Although we do not expect these relations to hold exactly, since they do not take into 
account the boundaries in z, it can be seen that they are reasonably well satisfied; £ x is 
approximately linear in Q (the discontinuity in the derivative of £ x simply reflects the fact 
that £ x is odd in Q whereas £ x is an rms quantity), whilst £ y (which is even in Q) shows 
some variation about a significant mean value. This reflects what is expected from the 
discussion of symmetries in §(2j where we explained how we would expect a change in the 
sign of Q to lead to a change in the sign of £ x but not of £ y . 



In the traditional formulation of mean field electrodynamics (see iMoffattl 119781 ). 
regeneration of the mean field relies on a lack of reflectional symmetry in the flow, usually 
characterised by the flow possessing net helicity. Field regeneration is then ascribed to 
the 'a-effect', which is given by the symmetric part of the Oj,- tensor in expression Q; 
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the antisymmetric part of ajjj is interpreted as a mean pumping velocity for the magnetic 
field. As discussed above, the decomposition (E]) is not meaningful for the problem we are 
considering. However it is still meaningful to examine which parts of the emf depend on 
reflectionally symmetric parts of the physics. This is most readily achieved by changing 
the sign of Q - the elements of £ that are reflectionally symmetric will change sign with Q 
while the non- reflectionally symmetric parts will not. In the absence of any rotation, £ x is 
zero while £ y is at its largest in magnitude. This leads us to associate £ y with magnetic 
pumping and turbulent diffusion. 
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Fig. 3. — Variation of £ x , £ y and s r with Q. All other parameters are as in Fig. [TJ 



Figure |H displays the dependence of £ x , £ y and s r on colatitude 9; note that the 
instability is more vigorous at the poles. Using only simple symmetry arguments, 
equation (1281) suggests that £ x oc cos 9 and £ y = constant. The agreement actually is quite 
good at high latitudes, but rather less so for mid-to-low latitudes. 

Figure |5] shows the variation of £ with the basic state magnetic field gradient (. It can 
be seen that the growth rate increases approximately linearly with (, confirming the idea 
that a larger field gradient leads to a more vigorous instability. The magnitude of both 
£ x and £ y increases with £ for small values, which is to be expected as the instability sets 
in and becomes stronger. In general, the significance of £ x reflects a balance betwe en the 
vigour of the instability and the influence of rotation, as illustrated by the picture of [Parker 
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Fig. 4. — Dependence of S x , S y and s r on 9. All other parameters are as in Fig. [TJ 



( 119551 ) of rising, twisting flux loops (termed 'cyclonic events' by Parker). If the instability 
is weak then S x is also weak; conversely if the instability is strong, the influence of rotation 
is diminished and £ x is reduced. This feature can be seen clearly in the plot of £ x versus (. 
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Fig. 5. — Dependence of £ x , £ y and s r on (. All other parameters are as in Fig. [TJ 



Figure E] shows the variation of £ x , S y and s r with A, which corresponds to a change in 
the strength of the basic state magnetic field at the top of the layer (since p(z = 0) is kept 
constant). For low values of B , S x increases — this can be thought of as an increase in the 
field strength allowing the instability to become stronger and therefore generate more emf. 
As the strength of the field increases further, the growth rate continues to increase; the 
influence of rotation is thus diminished and, just as in Figure |5j S x decreases in magnitude. 
There is a sharp increase in S y over low Alfven speeds as the instability sets in, a slight 
decrease as the Alfven speed is further increased and then a sharp drop at high Alfven 
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5. Sheared Magnetic Field 

Since for many astrophysical dynamos it is believed that the toroidal field, generated 
by the differential rotation (the 'w-effect' in the language of mean field electrodynamics), 
dominates the poloidal component, then it makes sense initially to consider the nature of 
the emf that results from the instability of a purely toroidal field; this is the approach 
adopted in §§[3],IU The dynamo resulting from considering only the w-effect and the mean 
emf resulting from a unidirectional field may be considered as the analogue of the classical 
ctcu-dynamo. However, in reality, the generation of a mean y-component of the field will 
influence the instability to some degree, and it is therefore of interest to extend the analysis 
of §4 so as to consider the nature of the emf arising from the instability of a sheared 
horizontal magnetic field. In this section, therefore, we consider a basic state magnetic field 
of the form B = B x e x + B y e y , where in non-dimensional form B x = (1 + (z) as before, and 
B y = A(l + £z). The extended problem thus has two additional parameters: £, the gradient 
of the ^-component of the basic state magnetic field, and A, the ratio of the ^-component of 
the basic state magnetic field to the x-component at the top of the layer. 

The influence on the dynamics of the instability of a weak y-component of the magnetic 
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field is of interest in itself, with applications to the morphology of the field that may emerge 
from the solar tachocline. The inclusion of B y leads not only to changes in the stratification, 
but also to an additional restraint on the instability, since any instability must 'unwind' 
the field. We propose to discuss the details of the instabi lity in a separate pap er; here we 



Cattaneo et al. 



(11989h have also 



shall concentrate on its influence on the emfs generated, 
considered magnetic buoyancy instability of a sheared magnetic field, but they restricted 
attention to the non-rotating case; in order to determine a preferred ('toroidal') direction 
they made the restriction to x-independent perturbations (k — 0). In our system, the 
presence of rotation distinguishes the x- and y-directions and we may therefore consider 
fully three-dimensional disturbances, without ambiguity. 

Figures [7] - [9] illustrate the dependence of £ on various parameters, with emfs weighted 
using the same prescription as in §|U Unlike in §HJ however, the two fastest-growing modes 
no longer necessarily have oppositely-signed x-direction wave numbers. For a unidirectional 
magnetic field, only the angular rotation O breaks the symmetry between the (k, I) and 
(—k, I) modes. The introduction of a poloidal ingredient means that there is an additional 
symmetry breaking, and so there is even less reason to expect the fastest growing modes to 
have oppositely-signed k. The weighting procedure detailed in §|2]is still applied so that the 
results can be compared directly in the case of B y — Y with those of §HJ 

Figure [7] displays the variation of £ x , S y and s r with Q. The growth rate and both 
components of the emf are diminished as \Q\ is increased. The growth rates of the two 
fastest-growing modes are very similar, as can be seen from the fact that the solid and 
dotted lines in the figure are indistinguishable. Unlike in Figure |3j the addition of the 
poloidal component B y means that S x is no longer antisymmetric with respect to Q, and, 
although from the graph it is not discernible, S y is no longer symmetric in Q. With non-zero 
By, the non-rotating case no longer represents a significant case of symmetry; when B y = 0, 
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the case of f2 = represents a situation in which the poloidal and toroidal directions cannot 
be distinguished — the introduction of the poloidal component removes this property. Now, 
£ x has both reflectionally symmetric and non-reflectionally symmetric parts, as does £ y . 
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Fig. 7. — Dependence of £ x , £ y and s r on Q in the case with a sheared basic state magnetic 
field. Parameter values are 9 = 75°, A 2 = 0.1, ( = l,v = rj = 5x 1CT 5 , x = 2, A = 0.1 and 

Figure El shows the variation of £ x , £ y and s r with A (the ratio B y /B x at the top of 
the layer). As we have chosen a case with £ > 0, both the x- and y-components of B are, 
individually, potentially destabilising (field strength increasing with depth); increasing A 
therefore increases an already destabilising component of the field. As the strength of the 
poloidal field increases, the growth rate increases, as does the emf in the x-direction, while 
the emf in the y-direction slightly decreases in strength. The increase in the value of £ x can 
be viewed as the instability becoming stronger as A is increased. As expected, the difference 
between the two fastest-growing modes increases with B y , as is evident from the separation 
of the two lines in Figure [B^c). 

Figure |9] shows the dependence of £ x , £ y and s r on £, the gradient of B y , for £ = 1. 
Both positive and negative values of £ are considered — positive values of £ correspond to a 
basic state field in which both the x- and y-components can be thought of as destabilising, 
while negative values of £ correspond to a stabilising .E^-component. When £ < 0, both 
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Fig. 8. — Dependence of £ x , £ y and s r on A at the top of the layer. Parameters as in Fig. [7J 
with O = 0.1. In the third panel, the dashed and solid lines represent the growth rates of 
the two fastest-growing modes, which are then used in the weighting process. 

the emf and the growth rate increase as £ becomes less negative (i.e., as the poloidal 
field increases less with height), and then continue to increase with £ for low values of £. 
Eventually the emfs and growth rates all begin to decrease as £ increases further. 
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Fig. 9. — Dependence of S x , £ y and s r on £. Parameters as in Figs. [7J and |HJ 



In their study of the magnetic buoyancy instabilit y of a sheared magnet ic field, with 
the crucial restriction to ^-independent perturbations, ICattaneo et al.l ( 119891 ) found that 
the nature of the instability was strongly influenced by the location of the height at which 
By = 0, where the field appears locally untwisted to x- independent perturbations (the 
resonant surface). It is therefore of interest to examine this feature in our system, although 
it is important to note that here we are considering fully three-dimensional instabilities, for 
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which the notion of a unique untwisted direction is not clearly defined. We have examined 
the effect of varying the depth at which B y = (henceforth called z ), whilst keeping 
the magnetic energy of the ^-component of the field constant. The ratio of the magnetic 
energies of the y- and x-components of the basic state magnetic field is 3/7. 

Figures [TU], [TT] and [T2] show £ x and £ y as functions of depth for three values of z : 
Zq — 0.1 {By = near the top of the layer), z = 0.5 (B y = in the middle of the layer) and 
zq = 0.9 (By = near the bottom of the layer). Results with z < or z > 1 are extremely 
similar to those with zq near the top or bottom of the layer respectively, and so are not 
displayed here. The quantities displayed are £ x and £ y , calculated using the eigenfunctions 
of the single most unstable mode, unweighted so as to allow direct comparison with Figured] 
(which is an example of an unweighted emf resulting from a unidirectional basic state field). 

It can be seen that the components of the emfs in each case take very small values 
throughout much of the layer; the depths at which £ x and £ y are significant though depend 
on Zq. For all values of zq, the emfs assume a similar shape to those seen in Figure [TOl 
Comparing the three plots of £ x in Figures fTOl - fT2l we can see that as z$ is increased, the 
first peak (with £ x > 0) decreases in height, while the second (with £ x < 0) is amplified. The 
variation of the distribution of the emfs within the layer is of interest. For a unidirectional 
field (Figure [T|), £ x and £ y are essentially non-zero only close to the top of the layer, in 
a similar manner to Figure [10J In cases where z lies either far above or far below the 
layer, the untwisted surface plays no role in the instability mechanism and so we expect a 
modified version of the unidirectional case. However, when zq is located within the layer, 
then the upper region of the layer has \B y \ increasing with height, and so has a stabilising 
influence, while the lower region has \B y \ decreasing with height, a destabilising influence. 
This suggests that the eigenfunctions should, in general, be dominant in the region z > zq, 
i.e. in the least stable region, as can be seen in Figure [TTJ When z Q is close to the bottom 
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boundary, this is not possible, owing to the boundary conditions, and so the eigenfunctions 
are non-zero towards the middle of the layer, as seen in Figure | 
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Fig. 10. — S x and S y for the the most unstable mode in the case z Q = 0.1. The wave numbers 
are k = 0.18987 and I = 7.4219, and the eigenvalue is s = 0.050992 - 0.043032i 
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Fig. 11. — £ x and E y for the most unstable mode in the case z Q = 0.5. The wave numbers 
are k = -0.35297 and I = 7.4875, and the eigenvalue is s = 0.02385 - 0.025888i. 



6. Conclusions and Discussion 

We have examined the nature of the mean electromotive force that arises from 
rotationally influenced magnetic buoyancy instability of a stratified magnetic field. In a 
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Fig. 12. — S x and 8 y for the most unstable mode in the case z = 0.9. The wave numbers 
are k = 1.6178 and I = 6.1094, and the eigenvalue is s = 0.021667 - 0.023000i. 

standard kinematic mean field dynamo the magnetic field is amplified from an infinitesimal 
strength via an a-effect that results from consideration solely of the velocity field; an 
instability-driven dynamo, on the other hand, is intrinsically dynamic, in that the motions 
responsible for driving the emf arise only from the instability of the mean magnetic field 
itself. Although instability-driven emfs have been studied in the context of magnetic 
confinement devices, they have received relatively little attention in the astrophysical 
literature. However, magnetically-driven instabilities arise in a number of astrophysical 
contexts and, for this reason, together with possible difficulties associated with the standard 
mean field approach, we believe it is of some interest to investigate this idea in some detail. 

One of the most important points to note is that for an instability-driven emf, the 
standard decomposition fl3]) is no longer meaningful and that it is the total emf that is of 
interest here. It is though still possible to decompose the emf into components that do and 
do not depend on a lack of reflectional symmetry in the flow. The dependence of the emf on 
the mean field and its gradient can be extremely complicated, as shown by expression (137]) . 
which considers the relatively simple case of an ideal MHD analysis with v = rj = 0. For 
the more general treatment of the instability, discussed in §§HJ[5J the stability problem and 
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the subsequent evaluation of the emf have to be performed numerically. In order to treat 
the competition between modes of equal growth rates, but different emfs, it is necessary 
to introduce some weighting procedure, such as that described by expression fTSBl) . We can 
then obtain both the latitudinal and depth dependence of the emfs, as well as exploring the 
dependence on parameters such as rotation. 

We have shown that a net emf arises naturally from the preferred modes of magnetic 
buoyancy instabilities in a rotating system. Our results also make clear that the relationship 
between the effectiveness of the emf and the vigour of the instability is not straightforward 
- this reinforces the important message that an instability per se cannot guarantee a mean 
emf and hence a closing of the dynamo loop. 

Our work described here has addressed only the nature of the mean emf that arises 
from magnetic buoyancy instability. The crucial next stage is to investigate the viability of 
the wider dynamo process, by incorporating such emfs into a model of the evolution of the 
mean field. We propose first to consider a highly simplified model, with the emfs taking 
prescribed functional forms of the mean field and its derivatives, based on the model in 
§|3l before considering an extended model in which the instability of the mean field and 
its evolution under the resulting emf are cal culated seff-consisten tly. This latter approach 



has some similarities with that adopted by 



Turner et al 



(120091 ) . who considered the 



hydrodynamic problem of the diffusion of the mean vorticity profile, through consideration 
of the products of linear fluctuations; a crucial difference though is that in their problem 
the response was wavelike, not the result of an instability. 

Even in the broadest sense, there is currently no consensus of how t he solar dynamo 



opera tes. Three possible scenarios have been advanced (see, for example, 



Tobias fc Weiss 



20071 ): a distributed dynamo, in which field regeneration takes place throughout the 



convection zone; an interface dynamo, in which all the action is concentrated at the base 
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of the convection zone and in the tachocline; and a flux transport dynamo, which might 
be envisaged as having an cu-effect arising from the velocity shear in the tachocline and 
an ct-effect located close to the solar surface, with the two linked by some large-scale flow. 
There are certainly problems with all three types of dynamo. A distributed dynamo relies 
on a cohere nt cc-effect throughout the convection zone; this is no t supported by th eoretical 



arguments (jCattaneo fc Hughes 



20091 ) or numerical simulations ( jBrun et al. 



20041 ). A flux 



transport dynamo relies crucially on the migration of surface magnetic features and a 
large-scale flow to return the field to the base of the convection zone, essentially ignoring 
all the turbulent dynamics of the convection zone. An interface dynamo is more appealing 
in that the different aspects of the dynamo process are not widely spatially separated. 
Toroidal field would be amplified from the poloidal component by the differential rotation in 
the tachocline; the difficulty though is in closing the dynamo cycle. One possibility is from 
helical overshooting convection, though this would then still rely on the traditional a-effect 
working at high magnetic Reynolds numbers. An emf resulting from magnetic buoyancy 
instability provides a natural possible alternative solution to this problem. Clearly, further, 
more involved, investigations are needed in order to examine the viability of this idea. 
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